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ABSTRACT 

Misfolding and associated loss of function are 
common problems in constructing fusion RNA 
complexes due to changes in energy landscape 
and the nearest-neighbor principle. Here we report 
the incorporation and application of the pRNA-3WJ 
motif of the phi29 DNA packaging motor into fusion 
RNA with controllable and predictable folding. 
The motif included three discontinuous ~18 nucleo- 
tide (nt) fragments, displayed a distinct low folding 
energy (Shu D ef a/., Nature Nanotechnology, 2011, 
6:658-667), and folded spontaneously into a leading 
core that enabled the correct folding of other 
functionalities fused to the RNA complex. Three 
individual fragments dispersed at any location 
within the sequence allowed the other RNA func- 
tional modules to fold into their original structures 
with authentic functions, as tested by Hepatitis B 
virus ribozyme, siRNA, and aptamers for malachite 
green (MG), spinach, and streptavidin (STV). Only 
nine complementary nucleotides were present for 
any two of the three ~18-nt fragments, but the 
three 9 bp branches were so powerful that they dis- 
rupted other double strands with more than 15 bp 
within the fusion RNA. This system enabled the pro- 
duction of fusion complexes harboring multiple RNA 
functionalities with correct folding for potential ap- 
plications in biotechnology, nanomedicine and 
nanotechnology. We also applied this system to in- 
vestigate the principles governing the folding of 
RNA in vivo and in vitro. Temporal production of 
RNA sequences during in vivo transcription caused 
RNA to fold into different conformations that could 
not be predicted with routine principles derived 
from in vitro studies. 



INTRODUCTION 

RNA nanotechnology involves programmable and 
addressable designs of RNA 3D nanoparticles by fusing 
or grafting individual RNA structural moieties or 
functionalities to one another (1,2). Many RNA folding 
programs are available (3-5) for RNA secondary structure 
prediction. However, manipulating tertiary folding of 
the reconstructed RNA is very challenging due to changes 
in the nearest-neighbor and the redistribution of energy 
landscapes. Development of systems to fold RNA 
functionalities in a controllable and predictable manner 
is therefore very desirable. It has been reported that RNA 
junction motifs, such as the three-way junction (3WJ) (6-9), 
four-way junction (4WJ) (10), and five-way junction (5WJ) 
(11) can form stable conformations composed of discon- 
tinuous RNA fragments. For example, the stable properties 
of the 3WJ structural motif derived from pRNA of bac- 
teriophage phi29 DNA packaging motor has been utilized 
to generate functional RNA nanoparticles (6-8,12). Other 
examples of in silico generated RNA nanostructures based 
on RNA structural motifs include RNAI/II inverse kissing 
loops (13), kink- turn structural motifs (14,15), receptor- 
loop/loop-receptor interacting motifs (16,17), and phi29 
pRNA hand-in- hand and foot-to-foot interactions (1,12, 
18-22). In addition, tightly folded RNA motifs have been 
reported in vivo that served as vectors to carry exogenous 
functional RNA modules (23). 

More and more non-coding RNA molecules have been 
discovered to play essential roles in the regulation of a 
variety of diverse biological functions in cells (24-27). 
Often, the function of an RNA molecule is dependent on 
its tertiary structure, as well as information encoded in its 
canonical or non-canonical base pairing (28-30). Thus, 
prior to in vitro or in vivo design of artificial RNA 
nanoparticles, one needs to consider if the desired fusion 
RNA will fold into the predicted functional structure and 
retain its correct functionality. An example is the RNA 
paranemic motif coupled with malachite green (MG) 
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(a) 3WJ motif (b) 




Figure 1. Experimental design for the construction of fusion RNA 
complexes harboring multiple functionalities driven by pRNA-3WJ 
core, (a) Secondary structure of the pRNA from the bacteriophage 
phi29 DNA packaging motor. The central 3WJ motif is boxed, 
(b) Approach for the construction of fusion RNA complexes. In the 
figure, a 3WJ , b 3WJ , and c 3WJ represent the three strands of the pRNA- 
3WJ complex; 1—4 represent RNA functional modules. 

binding RNA aptamer (31) to confer rigidity and enhance 
fluorescence emission (6,7,8,32). 

The structural properties encoded in intracellular func- 
tional RNAs are extremely important in RNA nanotech- 
nology. Several methods have been applied to determine 
structural properties of small and large RNAs (33-37), 
including physical and chemical approaches to probe 
RNA/ribonucleoprotein structures in vitro (38). 
Additionally, several chemical reagents that are sensitive 
to secondary and/or tertiary structures have been used for 
elucidating RNA structure in vivo (34). Dimethyl sulfate 
(DMS) has been extensively applied for probing RNA 
structure in a variety of organisms ranging from bacteria 
to eukaryotes (39,40). Lead-(II)-acetate has also been used 
to probe RNA structures in bacteria (41); this ion easily 
enters bacterial cells and primarily induces specific cleav- 
ages at positions of tight metal ion binding. Another 
chemical reagent is the hydroxyl radical (OH), which has 
been widely used in nucleic acid footprinting (42). This 
technique can be applied to probe RNA tertiary structure 
and its intermolecular interactions. 

Recently, we discovered that the pRNA-3WJ motif 
(Figure la) exhibits unusually robust properties (6,8,43) 
and can serve as a core to drive the folding of other 
RNA molecules fused to the complex. The pRNA-3WJ 
domain can be assembled from three pieces of RNA oligo- 
nucleotides resulting in a structure with: (i) an unusually 
high thermodynamic stability; (ii) highly efficient complex 
assembly even in the absence of metal salts; (iii) resistance 
to denaturation even in the presence of 8 M urea; (iv) the 
ability to stay intact without dissociating at ultra-low con- 
centrations; and (v) coordination of two divalent metal 
ions (6,43). Herein, we demonstrate that this robust 
pRNA-3WJ motif can provide a leading core and drive 
the correct folding of other functionalities fused into the 
RNA complex with controllable and predictable conse- 
quences. This system will enable the production of 
fusion RNA with incorporated functionalities folded 
into their authentic structure and function for 
applications in biotechnology and nanotechnology. 



previously (44). The DNA templates and primers were 
synthesized chemically by IDT (Iowa). RNAs were used 
directly for MG fluorescence assays, or purified by 8M 
urea 8% polyacrylamide gel electrophoresis (PAGE). 

Folate binding assay 

Human nasopharyngeal carcinoma KB cells [American 
Type Culture Collection (ATCC)] were grown on 
glass coverslides in folate-free medium overnight. Cy3- 
labeled pRNA-3WJ-FA-siRNA(survivin)-Ribozyme was 
incubated with the cells at 37°C for 2h. After washing 
with phosphate buffered saline (PBS), the cells were 
fixed by 4% paraformaldehyde and stained by Alexa 
Fluor® 488 phalloidin (Invitrogen) for cytoskeleton and 
TO-PRO®-3 iodide (642/661) (Invitrogen) for nucleus. 
The cells were then assayed for binding and cell entry by 
Zeiss LSM 510 laser scanning confocal microscope. 

Assay for the silencing of genes in cancer cell model 

KB cells were transfected with 25nM of individual 
fusion pRNA-3WJ complex using Lipofectamine 2000 
Transfection Reagent (Invitrogen). After 48 h, cells were 
collected and target gene silencing effects were assessed by 
western blot assay, as previously described (6). Briefly, 
cells were lysed by Radio-Immunoprecipitation Assay 
(RIPA) lysis buffer (Sigma) and the cell total protein 
was extracted for the assay. Equal amounts of proteins 
were then loaded onto 15% sodium dodecyl sulfate 
(SDS)-PAGE and electrophoretically transferred to 
Immun-Blot PVDF membranes (Bio-rad). The 
membrane was probed with survivin antibody (R&D) 
(1:4000 diluted) and (3-actin antibody (Sigma) (1:5000 
diluted) overnight, followed by 1:10000 anti-rabbit sec- 
ondary antibody conjugated with horseradish peroxidase 
(Millipore) for 1 h. Membranes were blotted by ECL kits 
(Millipore) and exposed to film for autoradiography. 

Binding assay for the Streptavidin-binding aptamer 

Tritium-labeled [ 3 H] RNA nanostructures (pRNA-3WJ- 
MG-HBV-STV) were preassembled in binding buffer 
(PBS with lOmM Mg 2+ ) before incubation with 
streptavidin (STV) agarose resin (Thermo Scientific). A 
total of 50 ul aliquot of STV resin was equilibrated at 
room temperature (RT) following washing with binding 
buffer. In all, ~5ug RNA samples were added to each 
tube and incubated with the resin for 1 h at RT. After 
incubation, the resin was spun at 500 x g for 1 min, and 
the supernatant (pass through) was removed. Then, 50 ul 
binding buffer was added to the mixture and incubated for 
15 min to wash the resin several times. RNA was eluted by 
5mM biotin solution and samples were analyzed on a 
1900 TR Liquid Scintillation Counter (Packard). 



MATERIALS AND METHODS 

In vitro synthesis and purification of RNAs 

All pRNA-3WJ constructs were prepared by in vitro tran- 
scription using T7 RNA polymerase from DNA templates 
amplified by polymerase chain reaction (PCR), as described 



Hepatitis B virus ribozyme activity assay 

Hepatitis B virus (HBV) substrate was labeled with Cy3 
(Mirus Bio LLC) and incubated with RNA nanoparticle 
at 37°C for 1 hr in Tris buffer (50 mM Tris-HCl, 
pH = 7.5; 20 mM MgCl 2 , 20 mM NaCl). pRNA-HBV 
ribozyme was used as a positive control (45). After 
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incubation, the products were run in 8M Urea, 10% 
PAGE gel for fluorescence imaging. 

MG aptamer fluorescence assay 

Gel-purified RNAs were mixed with 2 uM Malachite 
Green, tremethylmetane (MG) in binding buffer contain- 
ing lOOmM KC1, 5mM MgCl 2 and lOmM HEPES (pH 
7.4) and incubated at RT for 30min. The refolded RNA 
samples were heating to 95°C for 5 min followed by slow 
cooling to 37°C prior to MG staining. Fluorescence was 
measured using a fluorospectrometer (Horiba Jobin Yvon; 
SPEX Fluolog-3), excited at 615 nm, scanning from 625 to 
800 nm for emission (7,32,46). 

Spinach fluorescence measurements 

Assembled pRNA-3WJ nanoparticles (0.2 uM) fused with 
Spinach aptamer in TMS buffer were mixed with DFHBI 
dye (2 uM) and incubated at room temperature for 30 min. 
Fluorescence was measured using a fluorospectrometer 
(Horiba Jobin Yvon), excited at 450 nm and scanning 
from 565 to 750 nm for emission. 

Vector construction for in vivo expression 

A sequence of cis-ribozyme was fused onto the 3'-end of 
DNA templates of corresponding RNA for terminal pro- 
cessing (45,46). The 5' + 0 (Reference), 5' + 6, 5' + 12, and 
5' + 15 were inserted between Bglll/Ndel sites of expression 
vector pET-3b. Bglll cleavage removed the original T7 
promoter in the vector to prevent any undesired sequence 
in the 5'-end of the RNA product. The insertion fragment 
contained the T7 promoter at the 5'-end. The clone was 
transformed into Escherichia coli strain DH5a and recom- 
binant plasmids were verified by sequencing (GENEWIZ). 

RNA in vivo expression and purification 

Escherichia coli strain BL21 star (DE3; Invitrogen) was 
transformed by the recombinant plasmids, and the 
colony was inoculated by 5 ml LB medium containing 
lOOug/ml ampicillin, grown at 37°C, and shaken at 
250 rpm until A 60 o nm reached 0.5. IPTG solution of 
50 ul (1 M) was added to 5 ml of cell culture. Cells were 
allowed to grow for 1.5 h, and then were pelleted and re- 
suspended in 250 ul of 10 mM magnesium acetate, 1 mM 
Tris-HCl, pH 7.4 (47). The total soluble RNAs were ex- 
tracted using 500 ul of water saturated phenol (pH 4.5; 
Fisher). The aqueous phase was ethanol precipitated 
then dissolved in 50 ul of 0.05% DEPC treated water. 

Gel electrophoresis of in vivo prepared RNAs 

In vivo expressed 5' + 0 (reference), 5' + 6,5' + 12, and 5' + 15 
RNAs were analyzed by 8% PAGE gel containing 8M of 
urea by loading 5 ul of each sample and run in lx TBE 
(89 mM Tris base, 89 mM Boric acid, 2 mM ethylenediami- 
netetraacetic acid) at 120 V, at RT for 1 h. The 8% native 
PAGE gel was run in lx TBM (89 mM Tris base, 200 mM 
Boric acid, 5mM MgCl 2 , pH 8.0) at 70 V, at 4°C for 3h. 
After electrophoresis, the gels were stained with 10 uM MG 
in 100 mM KC1, 5mM MgCl 2 , and lOmM HEPES (pH 7.4) 
for 15 min at RT. The MG fluorescence image was acquired 



by the Cy5 channel (635 nm excitation/670 nm emission) 
using a Typhoon scanner (48). The gels were then stained 
with EB and scanned in EB channel (532 nm excitation/ 
580 nm emission). 

RESULTS 

pRNA-3WJ motif drives the folding of different RNA 
functionalities in fusion RNA 

Figure lb shows the design principles for constructing 
fusion RNA with multiple functionalities. The three 
pRNA-3WJ fragments (a 3 \yj, b 3W j, and c 3W j) are 
dispersed within the complex and serve as a driving 
force for the folding of other RNA molecules. Several 
RNA functionalities were incorporated and tested 
in vitro including: (i) siRNA targeting survivin gene 
(49,50); (ii) MG (Malachite Green dye, triphenylmethane) 
binding RNA aptamer (51); (hi) Spinach RNA aptamer 
(52); (iv) STV binding RNA aptamer (53); and (v) HBV 
ribozyme (45). These RNA moieties are unique in their 
sequences and folding properties, and exhibit different 
Gibbs free energies (AG). When the functional therapeutic 
RNA moieties, i.e. siRNA, ribozymes or receptor-binding 
aptamers, were fused to any of the three branches of the 
pRNA-3WJ, the 3WJ and the incorporated RNA modules 
folded independently, as validated by AFM image (Figure 
2a), and retained their authentic functions, as 
demonstrated by functional assays (Figure 2). The 
location of the three 3WJ fragments in the fusion RNA 
was not critical; but the order of appearance from 5'- a 3W j, 
to b 3WJ , to c 3WJ -3' (Figure lb) was the essential factor. 

Assessment of cell binding and gene silencing using fusion 
pRNA-3WJ-FA-siRNA (survivin)-Ribozyme nanoparticles 

A pRNA-3WJ complex was coupled with folate (FA) 
ligand, survivin siRNA, and ribozyme, denoted as 
pRNA-3WJ-Folate-siRNA(survivin)-Ribozyme nanopar- 
ticle (Figure 2a). FA serves as a cancer cell delivery 
agent via receptor mediated endocytosis (6,8,54). Cy3- 
labeled pRNA-3WJ complex was incubated with FA- 
receptor positive KB cells and tested for cell binding effi- 
ciency. Confocal microscopy indicated strong binding of 
RNA nanoparticles as demonstrated by excellent co-local- 
ization of cytoplasma (green) and RNA nanoparticles 
(red) (Figure 2b). 

Western blot was performed to assay the gene silencing 
effects of survivin siRNA in KB cells. After 48-h transfec- 
tion, reduced survivin protein expression level was 
observed compared to scramble controls (Figure 2c). 
P-actin was used as an endogenous control. We previously 
demonstrated that the mechanism of siRNA release from 
pRNA nanoparticles is by Dicer processing (21). 

Assessment of aptamer functionalities using fusion pRNA- 
3WJ-MG-Spinach nanoparticles 

MG binding RNA aptamers (51) and Spinach RNA 
aptamers (52) were fused to the pRNA-3WJ scaffold 
(denoted pRNA-3WJ-MG-Spinach nanoparticles) for 
structure and function verification (Figure 2d). Free 
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Assay for cell binding and targeted gene silencing 
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Figure 2. Functional assays of fusion RNA complexes harboring 
multiple functionalities, (a) Construction of RNA complex harboring 
folate (FA), survivin siRNA, and ribozyme (denoted, pRNA-3WJ-FA- 
siRNA (survivin)-Ribozyme) and corresponding AFM image, (b) 
Confocal images showing the targeting of FA-receptor positive KB cells 
by colocalization of cytoplasm (green) and RNA nanoparticles (red). Blue 
represents nuclei, (c) Western blot assay showing target knockdown of 
survivin protein by siRNA. P-actin serves as endogenous control, (d) 
Construction of RNA complex harboring malachite green (MG) (6- 
8,32,51) binding and Spinach aptamer (denoted, pRNA-3WJ-MG- 
Spinach) (52). (e) The 8% native PAGE verifying the fluorogenic 
properties of the two RNA aptamers. (f) Fluorescence spectra of MG 
aptamer (top) and Spinach aptamer (bottom) in solution, measured by 
fluorospectrometer. (g) Construction of RNA complex harboring MG 
aptamer, HBV ribozyme, and streptavidin aptamer (pRNA-3WJ-MG- 
HBV-STV). (h) The 8% native PAGE gel showing the catalytic activity 
of the incorporated HBV ribozyme. (i) MG aptamer activity verified by 
MG staining, (j) EB staining to verify the presence of RNA nanoparticles. 
(k) STV binding assay using STV affinity column (53). Bound RNA 
complexes were eluted by adding 5 mM biotin (outlined in red). 

MG and DFHBI (3,5-difluoro-4-hydroxybenzylidene 
imidazolinone; Spinach aptamer-binding dye) are not 
fluorescent by themselves, and only emit fluorescence 
after binding to the correctly folded RNA aptamers. 
Both MG binding aptamer and Spinach aptamers 
incorporated in the 3WJ motif retained their capacity to 
bind MG and DFHBI dyes, respectively, as demonstrated 
by their fluorescence emission in native PAGE (Figure 2e) 
and in solution (Figure 2f). 

Assessment of ribozyme catalytic activity and STV 
binding using fusion pRNA-3WJ-MG-HBV-STV 
nanoparticles 

Tetravalent RNA nanoparticles were constructed using 
pRNA-3WJ motif harboring MG binding aptamer 



(6,8,32), HBV ribozyme (45), and STV aptamer (53,48) 
(denoted pRNA-3WJ-MG-HBV-STV nanoparticles) 
(Figure 2g). HBV ribozyme (45) was able to cleave its 
135-nt RNA genome substrate into two fragments (60 nt 
and 75 nt) (Figure 2h). The results are comparable to 
optimized positive controls, thereby confirming that the 
fused ribozyme retained its authentic structure and 
function. MG aptamer also retained its capacity to bind 
MG dye, as revealed by fluorescence emissions in urea 
PAGE (Figure 2i) and verified by total RNA staining 
with EB (Figure 2j). 

For assessment of STV binding, [ 3 H]-UTP whole chain 
labeled pRNA-3WJ-MG-HBV-STV nanoparticles were 
incubated with STV resin in an affinity column. The 
column was washed, and the RNA nanoparticles were 
then eluted with biotin (Figure 2k), indicating correct 
folding of the incorporated STV aptamer. As a negative 
control, RNA nanoparticles without STV aptamer did not 
appear in the elution fractions. 



Design principles for evaluating the power of three 
pRNA-3WJ fragments to override other complementary 
sequences within fusion RNA constructs 

The energy landscape and nearest-neighbor principle are 
key factors governing RNA folding. Due to the free- 
energy minimization principle, complementary sequences 
within RNA also play a critical role in the formation of 
RNA 3D structures, e.g. pseudoknots, multi-way junc- 
tions, and hairpins. To apply the pRNA-3WJ for 
fabricating RNA nanostructures, the power of three 
pRNA-3WJ fragments to override other complementary 
sequences within fusion RNA complexes have to be tested. 
Accordingly, fusion RNA complexes were designed based 
on: (i) three fragments of pRNA-3WJ to direct the folding 
of the RNA; (ii) additional 5'-end interfering sequences 
complementary to the internal sequence of the 3WJ; and 
(iii) MG-aptamer to serve as a reporter for verification of 
authentic folding (7). If the RNA is degraded or 
misfolded, the MG fluorescence will disappear, and there- 
fore can be used as an indicator for assessing the folding of 
fusion RNA complexes. 

RNA constructs containing interfering sequences have a 
designated nomenclature. For instance, '5' + 12 means that 
12 overhanging nucleotides complementary to the MG- 
aptamer sequence were added at the 5'-end of the RNA 
complex. Similarly, "5'NM+12 (NM = non-matching) 
means that 12 nucleotides non-complementary to the 
MG-aptamer sequences were added at the 5'-end of the 
RNA complex. All RNA nanostructures were constructed 
to contain the pRNA-3WJ and an MG aptamer (Figure 
3a). Because the sequence at the 5'-end appeared earlier 
during transcription than the sequence at the 3'-end, 
single-stranded (ss) nucleotides of different lengths (6, 
12, and 15nt) were added to the 5'-end, complementary 
to the MG-aptamer sequence, in order to interfere with 
the original folding driven by the pRNA-3WJ (Figure 3b). 
Several controls were designed in which overhanging NM 
nucleotides of differing lengths complementary to the 
MG-aptamer sequence were added (Figure 3c). All these 
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(a) pRNA-3WJ-MG design 
scheme with modifications 



(b) pRNA-3WJ-MG with 5' end interfering sequences 



c) pRNA-3WJ-MG controls, with 5' end non-matching sequences 
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Figure 3. Design principles for evaluating the power of three pRNA-3WJ fragments to override other complementary sequences within fusion RNA 
constructs, (a) Schematic of the pRNA-3WJ structures showing location of the interfering nucleotides (left) and their complementary target (right). 
The 2D structures of the RNA complexes with 6, 12, or 15-nt interfering sequences (b) or scramble sequence (NM stands for non-matching) (c) at the 
5'-end (top panel) and corresponding fluorescence spectra (bottom panel). In (b) and (c), black and red curves represent fluorescence spectra after 
transcription and after refolding, respectively. 



RNA constructs contain discontinuous motifs, such as 
those observed in a pseudoknot (55-59). 

We then investigated whether the folding was affected 
by the discontinuous motifs prior to or after refolding 
subsequent to transcription. This is particularly important 
in the context of in vivo folding of RNA and for 
transcribing therapeutic RNA nanoparticles directly in 
the cell with unprocessed additional sequences at the 5'- 
or 3'-ends. Measurements of fluorescence signals of RNA 
complexes directly after in vitro transcription for all 
pRNA-3WJ constructs suggested folding of authentic 
MG-aptamer reporter system embedded within the 
pRNA-3WJ (Figure 3b and c, black curves). These 
results revealed that even with only nine complementary 
nucleotides within the 3WJ, the folding power of the three 
3WJ fragments could override more than 15 completely 
complementary nucleotides that could form double- 
stranded RNA within the RNA complex. 

In addition, a control experiment was performed in 
which the MG fluorescence intensity of the 5' + 12 con- 
structs was compared to the construct where the 3WJ 
motif was destroyed by deletion of the a 3WJ strand 
(Figure 4). The lower intensity of the control RNA 
complex suggests that the folding of the MG aptamer 
was indeed affected by the disruption of the 3WJ, 
thereby confirming that all 3WJ component strands need 
to be present to drive the folding of other RNA modules 
into their fully functional state. 

Comparison of pRNA-3WJ folding after transcription and 
after refolding 

To unravel the folding differences of in vitro RNA tran- 
scription products before and after refolding, pRNA- 
3WJ-MG emission spectra were measured. RNA tran- 
scription reactions were performed, as described previ- 
ously (60). The reactions were stopped by the addition 
of DNase I. Fluorescence spectra of the transcribed 
pRNA-3WJ-MG complex was then measured in the 
presence of MG fluorophore. The pRNA-3WJ-MG 



complex was then heated to 95° C and slowly cooled to 
37°C (refolding). The fluorescence spectra of the 
refolded pRNA-3WJ-MG complex was then compared 
with the transcribed complex. 

The fluorescence intensity was similar for all the control 
RNA samples, 5' + 0, 5'NM + 6, 5'NM+ 9, and 5'NM+ 12, 
regardless of the refolding step (Figure 3c). However, 
RNA constructs 5' + 12 and 5' + 15 with 12 and 15 over- 
hanging nucleotides, respectively, showed significant dif- 
ferences (Figure 3b). Upon increasing the ssRNA length 
at the 5'-end from 6 to 12 nt, as exemplified in molecule 
J' + 12, the fluorescence intensity after transcription was 
found to be lower than after refolding. Increasing the 
length of the ssRNA to 1 5-nt in RNA 5' + 15 resulted in 
a larger difference (Figure 3b). These results demonstrate 
that the complementary interfering sequences at the 5' -end 
partially disrupted the folding of pRNA-3WJ assembly 
during the transcription reaction. In contrast, if the 15- 
nt complementary interfering sequence is instead placed 
at the 3'-end, no differences in fluorescence intensity are 
observed (Figure 3C, complex 3' + 15), which is consistent 
with our interpretations. 

Ratio of fluorescence signals of pRNA-3WJ prior to and 
after refolding 

To assess the folding differences of the pRNA-3WJ prior 
to and after refolding, the integrated areas of individual 
fluorescence signals of RNA constructs were measured 
before and after refolding. The fluorescence ratio 'r' was 
obtained for each RNA sample by dividing the integrated 
area at 37°C (after transcription) by the integrated area 
after heating to 95°C (after refolding). These ratios were 
calculated from a minimum of four individual experiments 
(Table 1). The 'r' value reflects the impact of discontinu- 
ous motifs on the folding properties of pRNA-3WJ as a 
function of the length of the interfering sequence during 
transcription and after refolding (Figure 5). RNA con- 
structs with non-matching sequences to the MG-aptamer 
region, or with interfering sequences at the 3'-end, had 
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Figure 4. MG fluorescence intensity of the J' + 12 RNA construct com- 
pared to the construct where the 3WJ motif was destroyed by deletion 
of the a 3WJ strand, (a-b) The 2D structures of the 5' + 12 RNA com- 
plexes without (a) and with (b) deletion of a 3WJ strand, (c) Comparison 
of MG fluorescence spectra for RNA constructs in (a) and (b). 



Table 1. Ratio ('r') values of RNA nanostructures after transcription 
and after refolding 
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Figure 5. Correlation between the length of the interfering sequence 
and the ratio of fluorescence signals after transcription and after re- 
folding, (a) Descending trend of the ratio of fluorescence signals before 
and after refolding in accordance with the increasing length of the 5'- 
end interference sequences, (b) Control complexes with scrambled or 
non-matching sequences do not exhibit any change. 



similar 'r' values close to 1, indicating the effect was not 
dependent on nucleotide length (Figure 5b). However, in 
the presence of overhanging nucleotides that interfered 
with the MG-aptamer sequence there was a correlation: 
the longer the additional sequences at the 5'-end, the 
smaller the 'r' value, implying that there are fewer 
pRNA-3WJ-MG constructs that folded properly during 
temporal transcription (Figure 5a). The data leads to the 
conclusion that the RNA folding was temporally depend- 
ent. The 5'-end sequences that appeared earlier during 
transcription influenced the RNA folding more than the 
sequences that appeared later during transcription. 



In vivo assessment of the effect of temporal folding after 
transcription 

RNA constructs with interfering sequences of 6, 12, and 
1 5 nt, similar to the in vitro design of the 5' + 0, 5' + <5, 
5' + 12, and 5' + 15 RNA were expressed in E. coli cells. 
The principle of this design was based on the following 
ideas: (i) a thermodynamically stable and tightly folded 
pRNA-3WJ would drive the correct folding of the MG- 
aptamer; (ii) tightly folded RNA are relatively resistant to 
RNase degradation in the cell (most cellular RNases 
degrade ssRNA, but not dsRNA) and misfolded RNA 
would be trimmed or truncated by RNase, leading to a 
loss in fluorescence from the MG dye that only binds to 
the appropriately folded and intact MG aptamer; and 
(iii) the assay would be based on the comparison of the 
sequences that affect the folding of the core of the pRNA- 
3WJ-MG-aptamer. Thus, a cw-acting ribozyme sequence 
was placed at the 3'-end to cleave the undesired sequences 
arising from uncontrolled termination of transcription at 
the 3'-end (Figure 6a) (45,61). 

The RNA products from E. coli cells were first assayed 
by 8% denaturing PAGE (Figure 6b). As reported previ- 
ously, the pRNA-3WJ is resistant to 8M urea denatur- 
ation (6,8), and as such, denaturing gel can be used to 
detect the presence of the pRNA-3WJ-MG-aptamer 
complex. Upon staining with either ethidium bromide 
(EB) (Figure 6b, top) or MG dye (Figure 6b, bottom), 
the 5' + 0 RNA, serving as a reference molecule, and the 
5' + 6 RNA were detected on the gel as distinct bands 
around 200 bp, indicating the presence of both the 
pRNA-3WJ and the MG-aptamer in both constructs. 
The higher yield of 5' + 0 RNA compared to 5' + 6 con- 
struct indicated that the 6 complementary nucleotides at 
the 5'-end of 5' + 6 were able to disrupt the folding of the 
pRNA-3WJ-MG-aptamer complex only to some extent. 
However, MG-stained RNA was not detected in both 
5' + 12 and 5' + 15 RNA lanes, indicating the misfolding 
of those pRNA-3WJ-MG-aptamer complexes. Several 
degraded products were observed between 50 and 70 bp 
after EB staining (Figure 6b, top, lanes 5' + 12 and 5' + 75). 
This data agreed with the in vitro data, suggesting 
that RNA folding is temporally dependent during 
transcription, and that the 5'-end sequence appearing 
earlier during transcription influenced the RNA folding 
more than the sequence that appeared later during 
transcription. 

To verify the denaturing PAGE results, fluorescence 
intensities of the RNA constructs obtained in vivo were 
measured further in presence of MG dye (Figure 6c), in 
a manner similar to that of the in vitro experiments 
(Figure 3). The fluorescence intensities were detected for 
both 5' + 0 and 5' + 6 constructs (Figure 6c, black curves). 
But, the 5' + 0 RNA showed higher fluorescence value 
than the 5' + 6 construct with six interfering nucleotides. 
There was almost no fluorescence detected for constructs 
5' + 12 and 5' + 15. These data are consistent with MG- 
stained denaturing PAGE, suggesting that while six 
interfering nucleotides partially disrupted the folding of 
pRNA-3WJ-MG-aptamer, 12 and 15 nt complementary 
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sequences completely interfered with the formation of 
pRNA-3WJ-MG-aptamer structures in vivo. 

To refold the in vivo RNA products of pRNA-3WJ- 
MG-aptamer, the denaturation and refolding step was 
performed by heating the RNA solution to 95°C and 
slowly cooling to 37°C. Interestingly, the measured fluor- 
escence intensities were significantly reduced or were not 
detected at all (Figure 6c, red curves). These data are in 
sharp contrast to that obtained from the in vitro experi- 
ments. The plot of ratio versus nucleotide length from the 
in vitro RNA folding experiment resulted in a descending 
trend (Figure 5a), which was not observed for in vivo 
samples. The results suggest that RNA folding in vivo 
and in vitro are not always identical. Correctly folded 
5' + 12 or 5' + 15 RNA were significantly recovered after 
denaturation and refolding (Figure 3b). However, the 
misfolded in vivo 5' + 12 and J' + 15 RNA constructs, if 
any, could not be recovered upon refolding (Figure 6c). 

Comparison of the sequences with their complementary 
partners 

The folding properties of discontinuous 3WJ-pRNA were 
further analyzed by a prediction algorithm using free 
energy minimization parameters and mfold (4). The goal 
was to investigate how different lengths at the 5'-end 
contributed to the disruption of the 3WJ-pRNA con- 
structs, and whether there was a relationship between 
this influence and the predicted free energy of the added 
sequences on their partners. The following free energies 
for the duplexes of 5'- or 3'-end overhanging sequences 
and their complementary partners were calculated at 
10" 4 M RNA in the presence of 1 M NaCl: 

AG = - 6.1 kcal/mole for the 5' + 6 RNA 

AG = -19.4kcal/mol for the 5' + 12 RNA 

AG = -25.6kcal/mol for the 5' + 75 RNA 

AG = -26.6kcal/mol for the 3' + 15 RNA (control) 



The longest construct with a 15-nt complementary 
sequence had the highest stability with a AG of 
— 25.6kcal/mol. This indicates that the folding of the 
3WJ-pRNA complex 5' + 15 during in vitro transcription 
was strongly interfered. This theoretical data was consist- 
ent with the correlation observed by fluorescence measure- 
ments of 'r' values. The 3'-end interfering sequences 
resulted in the lowest free energy value, AG = 
— 26.6kcal/mol. This sequence did not interfere with the 
MG fluorescence signal since the transcription reaction 
occurs from the 5' to the 3'-end, and the 15-nt region 
was produced after the 3WJ complex had folded. 

DISSCUSION 

With the recent emergence of RNA nanotechnology (2), 
development of new systems to investigate the folding of 
RNA nanostructures has become very desirable. It has 
been reported that RNA junction motifs, such as 3WJs, 
4WJs, and 5WJs can form rigid structures composed of 
discontinued RNA fragments (6, 8,10,1 1). The application 
of robust RNA motifs as vectors to carry exogenous func- 
tional moieties has been reported (23,6,8,12,48). The 
novelty here is not the application of the RNA motif as 
vector, but the use of the unique pRNA-3WJ with favor- 
able thermodynamic attributes to drive the assembly of 
fusion RNA complexes (6,8). The energy landscape in 
RNA folding remains fundamentally important. When 
many motifs with similar levels of folding energies are 
fused into one large RNA molecule, new complementary 
sequences will appear and redistribution of energy land- 
scape following the nearest neighboring principle will 
occur. 

The pRNA-3WJ motif can be assembled from three 
pieces of RNA oligonucleotides with a T M slope close to 
90° (6,8). The resulting three-component complex was re- 
sistant to 8 M Urea denaturation and does not dissociate 
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at extremely low concentrations (6,8). The three compo- 
nents had a much higher affinity for favorable interactions 
compared to any of the two components, indicating co- 
operative simultaneous folding of the three helical stems. 
The incorporation of a pRNA-3WJ with distinct low 
folding energy will provide a leading motif to drive the 
folding of other functional motifs with weaker folding 
forces. Thus, functional modules with weaker folding 
forces will be able to fold independently into their 
original authentic structure when linked to the branches 
of the pRNA-3WJ. 

CONCLUSIONS 

The thermodynamically ultrastable pRNA-3WJ motif of 
phi29 motor pRNA can be used to drive the folding of 
fused RNA molecules with controllable and predictable 
consequences. When a large fusion RNA molecule is con- 
structed containing three individual 18-nt fragments 
dispersed at any location within the sequence, the three 
fragments drive the correct folding of other RNA 
functionalities or motifs into their original 3D structures 
with authentic functions. This approach offers a new 
system for future investigations into the important topic 
of RNA folding in vivo. 
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